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A major challenge in theoretical attempts to understand the microscopic 
mechanism of high Tc superconductivity is to quantitatively account for the su- 
perconducting condensation energy Microscopic model commonly used 
to describe the superconducting copper-oxides, the t-J model [0], is beleived 
to capture the essential physics of these materials: the interplay of electrons 
kinetic energy and their antiferromagnetic exchange interaction. Within the t-J 
model, the condensation energy can be related to the change in the dynamical 
spin structure factor between the superconducting and the normal states By 
analyzing previous experimental data, we show that the change associated with 
the resonant neutron scattering peak [p|- |ll|| found in Y Ba2Cu^0j can quantita- 



tively account for the condensation energy. We argue that this analysis suggests 
a microscopic mechanism for high superconductivity, where antiferromagnetic 
exchange energy is saved in the superconducting state through the coupling to 
a particle-particle resonance ||12[ -|16||. 

One reason for the absence of a generally accepted microscopic explanation of high 
superconductivity (HTSC) is the lack of a quantitatively precise definition of a "mechanism" . 
In the traditional BCS IQ superconductors, the phonon mechanism is quantitatively verified 
by the isotope substitution and single particle tunneling experiments. It is a widely held 
view that the HTSC involve strong electronic interaction, while phonon mediated interaction 
is not the major driving force for superconductivity. To date, many different "mechanisms" 
have been proposed, but very few quantitatively falsifiable predictions have been checked 
against experimental data. 
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A system undergoes a transition from the normal to the superconducting state because it 
can lower the total free energy. The condensation energy Ec-, defined as the energy difference 
between the normal state, extrapolated to zero temperature, {E^) and the superconducting 
state {Es)i can be directly measured through the thermodynamical critical field by the 
following relation [0: 

Ec = En — Es = -^Vq (1) 

OTT 

Here we choose Vo = ax6xctobe the volume of a unit cell (in the high Tc cuprates a = b 
with reasonably good accuracy) and define the energies inside one unit cell . Alternatively, 
the condensation energy can also be found by integrating the difference in specific heats in 
the (extrapolated) normal and superconducting states from T = to the superconducting 
transition temperature PJ^. 

On the other hand, the condensation energy can be computed in principle from a mi- 
croscopic Hamiltonian. Ideally, one would like to start from the full Hamiltonian involving 
electronic kinetic energy, the Coulomb interaction energy and the electron lattice coupling 
energy. In fact Chester |^ studied condensation energy of such a Hamiltonian and showed 
how various contributions to the condensation energy could in principle be determined ex- 
perimentally. Among many theories of HTSC, only the interlayer tunneling mechanism |^ 
made a quantitatively falsifiable prediction about the condensation energy. More recently, 
Leggett showed that the change in the interlayer Coulomb energy can be directly related to 
the frequency and momentum integral of the dynamic charge susceptibility, a quantity which 
can be measured experimentally 0. He argued that a microscopic theory of HTSC should 
quantitatively predict the frequency and the momentum range at which the change of the 
dynamical susceptibility occurs, and quantitatively account for the condensation energy. 

However, the energy scale involved in HTSC is much smaller than the Coulomb interac- 
tion. In order to understand the mechanism for HTSC, it is useful to start with an effective 
Hamiltonian which correctly describes the basic physics below the Coulomb energy scale. 
Recently, Scalapino and White |§] initiated the investigation of condensation energy in the 
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t-J model, which is defined by 



'H = -tY^clc,, + JY.SS, (2) 

Here is electron creation operator on site i with spin a, Sj is the electron spin operator, 
{ij) denotes a pair of near neighbor sites and the Hamiltonian acts on the space of no doubly 
occupied sites. Equation @) describes the Hamiltonian of one layer only. In a bi-layer system 
like Y Ba2Cu^0Q+x-, there are two layers per unit cell, with additional antiferromagnetic 
(AF) coupling J± between the layers. If a system described by this Hamiltonian undergoes 
a superconducting transition away from half-filling, the kinetic energy Et is expected to 
increase slightly, but the decrease in the exchange energy Ej may be significantly larger, so 
that a superconducting ground state is realized. 

In a tight binding model, the mean kinetic energy Et can be expressed as a frequency 
integral of the optical conductivity cr(a;) |]T7|-p!9[l. Then the contribution of Et to the conden- 



sation energy can be measured by comparing cr{uj) in the normal and the superconducting 
states. In the low frequency regime, the Drude peak in the normal state collapses into a 5 
function peak in the superconducting state, usually with conserved weight [^. Therefore, 
the change in the kinetic energy has to be determined from ct{uj) in the higher frequency 
range. 

Scalapino and White made a insightful observation that the change in the J term in 
equation (H) can also be directly expressed as a frequency and momentum integral of the 
dynamic spin structure factor S{q,oj). For bi-layer materials such as Y Ba2Cu-iOQj^x it is 
convenient to separate the even and odd parts of S{q., uo) with respect to bi-layer interchange 
within a unit cell 

S{q,uj) = S--{q,u)cos\^) + S"'\q,uj)sin\'^-^) (3) 

Here dcu is the distance between nearest-neighbor Cu-0 planes, so S'^^'^" and S°'^ describe 
the in-phase and out-of-phase spin fiuctuations in the two planes constituting one bi-layer. 
Then AEj = Ef — Ej may be obtained as ( per unit cell ) 



+ 



^ J (cosiq^a) + cos{qya)) + ^ 
^ J (cos(g^a) + cos{qya)) - ^ J 



} (4) 



From |2T|J2^ it is known that J± is much smaller than J, so the J± term in equation (^) may 
be safely neglected. Scalapino and White pointed out [§] that since S{q, u) can be measured 
directly in neutron scattering experiments, AEj can therefore in principle be measured. If 
AEj turns out to be much smaller than the experimentally measured condensation energy 
Ec, this gives a convincing way to rule out any mechanism of high Tc superconductivity based 
on the t — J model. On the other hand, if AEj turns out to be greater than Ec, one would 
have identified a important driving force for superconductivity, and can quantitatively test 
any theoretical mechanism of HTSC based on the t-J model, which should predict where in 
frequency and momentum space the change of S{q, uj) occurs. 

However, one should apply this line of reasoning with extreme care. The quantity 
SN{q,uj) in equation is not the normal state spin structure factor above Tc, but rather an 
extrapolated normal state quantity at T = 0. Experimentally, one has to carefully identify 
features in S{q, uj) which changes abruptly at Tc. Theoretically, one has to identify contribu- 
tions to S{q, Lu) which are absent in the normal state and only present in the superconducting 
state, so that the difference in equation (||) can be rigorously defined. Recently, a resonant 
neutron peak has been discovered in the family of YBa2Cu30Q+x superconductors [P-[TT[|. In 
the optimally doped Y Ba2Cu-iO'j superconductor with Tc = 91K, the resonance peak occurs 
in the spin scattering channel at energy ujq = 41 meV. It is centered around momentum 
Q = ^) and it occurs in the odd channel with respect to bi-layer interchange. Perhaps the 
most remarkable feature of the resonance peak is its onset at the superconducting transition 
temperature T^. Above Tc, there is no experimentally identifiable spectral weight in this 
frequency and momentum range. Therefore, this experiment is ideally suited for analyzing 
its contribution to the condensation energy. 

We constructed a theory |[T^ of the neutron resonance peak by first identifying a reso- 
nance in the spin-triplet particle-particle (p-p) channel defined by the following operator 



Tft = Y,{cosK - cosky)cl^Q^clk^ (5) 
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In a series of exact diagonalization studies on the Hubbard and t-J models, Meixner, Hanke 
and the two of us and Eder, Hanke and one of us (SCZ) demonstrated the existence 



of such a triplet p-p resonance for all doping range. A triplet p-p resonance can not be 
detected in the neutron scattering experiments on the normal state. However, below T^, 
the d-wave superconducting order parameter mixes a p-p resonance into the dynamical spin 
structure factor, leading to a sharply defined collective mode in the triplet particle-hole (p- 
h) channel (the so-called yr-resonance). The contribution from the p-p resonance to the p-h 
channel below Tc is given |12;1^ by 



d{hio) Siq,uj) = '^^^^^ (6) 
1 — n 

where (A^^) is the dimensionless d-wave superconducting order parameter (not the energy 
gap) and n is the density of electrons. This simple estimate of the spectral weight of 
the TT resonance at q = Q agrees reasonably well with the absolute intensity measured by 



Fong et al. [|rT|. On the other hand, the width of the vr resonance in momentum space 
can not be calculated reliably. Within the 5*0(5) theory [0], the tt resonance has been 
interpreted as a pseudo-Goldstone boson of the spontaneous S'0(5) symmetry breaking in 
the superconducting state. Its existence is therefore directly correlated with the existence 
of the superconducting long range order. 

Since the coupling to a p-p resonance is only possible in the superconducting state, it 
indeed contributes to the difference in S{q^uo) as required by equation (1). In this work we 
wish to point out that the logic leading to our explanation of the tt resonance below can 
be reversed to provide a microscopic mechanism of HTSC On one hand, the existence of 
the TT resonance is a unique property of the superconducting state, on the other hand, its 
existence around Q lowers Ej, as we see from equation (^). Therefore, the system undergoes 
a superconducting transition, so that a tt resonance mode can emerge, which in turn lowers 
Ej. In other words, no matter how antiferromagnetically correlated the normal state is, 
a superconducting state can always further lower the exchange energy by coupling a p-p 



resonance into the dynamic spin structure factor S{q,uj). Our mechanism for HTSC makes 
a quantitatively precise prediction which can be checked experimentally. Once we assume 
that the HTSC is dominantly driven by the emergence of a new collective mode, the tt 
resonance, this mechanism predicts that the condensation energy is determined by the zero 
temperature spectral weight of the vr resonance, and the energy saving occurs precisely in 
the frequency and momentum range of the vr resonance. 

Since our mechanism predicts an additional spectral weight around q = Q in the super- 
conducting state, a careful reader may wonder what happens to the spectral sum rule, which 
states that S{q,u!) integrated over both q and u (without the f{q) = cos^qxo) + cos{qyb) 
factor in equation (|D) should be a constant, independent of the nature of the ground state. 
Of course our additional spectral weight should be compensated by a depletion of spec- 
tral weights from other parts in the momentum space. However, since /(g) has a absolute 
minimum at g = Q, their contributions to (^ can not completely cancel the tt resonance 
contribution. In fact, depletion of the spectral weight near q = would lead to additional 
enhancement of the condensation energy. 

Since both the condensation energy and the spectral weight of the neutron resonance peak 
have been measured experimentally, it is straightforward to check this prediction. Within 
Landau-Ginzburg theory. He = „ %° where $o = ^ is the superconducting flux quantum, 
C, is the coherence length and A is the London penetration depth. In YBa2Cu30i, these two 
length scales are determined to be in the range of ^ = 12 — 20A and A = 1300 — 1500A. Using 
a = b = 3.85 A and c = 11.63 A this gives a condensation energy of Ec = 3.3 — 12 K per 
unit cell. On the other hand, Ec can also be measured directly in specific heat experiments. 
Loram et al. reported (see page 251 of [Q] and page 247 of 0) a value of Ec = 6 K per unit 
cell for YBa2Cu^07. Let us now turn to the experimental measurement of AEj given by 
equation (H). In practice, it is very hard to measure the change of S{q,u) throughout the 
frequency and momentum range. However, since our mechanism predicts that the major 
part of the condensation energy is saved in the range of ujq = 41 meV and a 2D momentum 
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transfer Q = (f,-^), we can analyze the neutron scattering data by assuming that the 
contributions of the other parts do not change significantly at Tc. Indeed, for frequencies 
above the resonance the neutron scattering intensity does not change as the system goes from 
normal to superconducting. For frequencies below and in the resonance range the normal 
state intensity Sn{q, i^) is below experimental sensitivity limit in both even and odd channels. 
Because of the sum rule discussed above we may not assume that it vanishes identically, 
so we make an assumption that in this frequency range the spectral weight Sn{q,^^) is 
spread uniformly in momentum space and therefore does not contribute to (|^). Below Tc, 
a neutron resonance peak emerges in the odd channel, and Sg'^'^{q,uj) has been measured in 
absolute units in the range of the resonance. The dimensionless quantity / d{Tiu)S'g^{Q^uj) 
is measured to be 0.52 at T = 10 K by Fong et al. (Theoretical estimate based on 
references [p!2|p!5| gives 0.32 for this value). The resonance has a Lorentzian profile centered 



aX q = Q with a width = 0.23A~^ |jTT|, so the 2D momentum integral can be easily 
estimated. 

3 / n \ 2 1 n S2 
AEj = - TT —K2D x-x^x2xJ = 0.016 J (7) 
2 V27r / 2 TT 

Taking J = 100 meV we obtain AEj = 18 K. The actual condensation energy should be 
smaller than AEj since the kinetic energy is expected to increase below the superconducting 
state. Taking this into account, we see that the spectral weight of the n resonance peak 
can quantitatively account for the condensation energy measured from and specific heat 
experiments, in reasonable agreement with the prediction based on our mechanism. 

It is worthwhile to compare our mechanism for HTSC with the spin fluctuation pairing 
mechanism [p^ -p5|. All these mechanisms are based on the AF exchange interaction inter- 
action Ej. Our analysis of the neutron data confirms that a large part of the condensation 
energy indeed arises from the AF exchange interaction 0]. However, our mechanism differs 
from spin fluctuation pairing mechanism, both conceptually and quantitatively. The idea of 
spin fluctuation pairing is directly borrowed from the phonon mediated pairing in the tradi- 
tional BCS superconductors. It requires sizable AF spin fluctuations in the normal state, in 



order to pair the fermi liquid like quasiparticles. In contrast, in our mechanism, the AF spin 
fluctuation does not play an important role in driving superconductivity, and is neglected in 
the normal state of optimally doped YBa2Cu^07 superconductor. In the superconducting 
state, the vr resonance can be viewed as a kind of spin fluctuation in a particular frequency 
and momentum range. From the point of view of the 5*0(5) theory [|T^, the superconducting 
state is obtained from the AF state by a rotation, it is therefore natural to expect it to have 
more AF correlation than the normal state. Therefore, the crucial conceptual difference 
between these two mechanisms lies in the fact that normal state AF spin fluctuations are 
not required in our mechanism. This conceptual difference has a direct quantitative conse- 
quence. From equation (|) we see that the presence of spin fluctuation spectral weight in the 
normal state contributes negatively to the condensation energy. It is not clear why the spin 
fluctuation model would predict a change in the frequency integrated weight of S{q, u>) near 
q = Q. In fact, most theories of the neutron resonance peak based on the spin fluctuation 
models assume a preexisting overdamped spin fluctuation mode in the normal state 

and only predict a reduction of damping below Tc, with essentially conserved weight. 

Our mechanism provides a natural explanation to the doping dependence of the conden- 
sation energy. We attribute the condensation energy to the difference in exchange energies 
in superconducting and normal states. More underdoped materials have considerably more 
AF correlations in the nomal state, as known from the neutron scattering p8[. Although 



in the superconducting state the weight of the resonance is enhanced the difference 
between the normal and the superconducting exchange energies becomes smaller. This is 
consistent with the results of Loram et al. who find a decrease in the condensation en- 



ergy with decreased doping in Yq,^C aQ,2Ba2CuzOQ+x |§3- Until this point we have been 
discussing the relation between the condensation energy and the resonant peaks in neutron 
scattering at zero temperature. It would be interesting to see whether this idea works at 
finite temperature as well. Without going into any details we would like to point at a strik- 
ing similarity between the temperature dependencies of H^. [0] and the resonance intensity 



|31|| for the underdoped cuprates. Not only both quantities scale similarly below Tc, but 
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they both have tails extending above Tc, presumablly arising from significant pairing fluc- 
tuations in the pseudogap regime of the cuprates. We postpone the quantitative analysis of 
these phenomena to future investigations. 

We would like to thank P. Dai, H. Fong, S. Heyden, B. Keimer, H. Mook, and D. Scalapino 
for useful discussions. 
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